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ABSTRACT 

Disc accretion to a rotating star with a non-dipole magnetic field is investigated 
for the first time in full three-dimensional (3D) magnetohydrodynamic (MHD) simula- 
tions. We investigated the cases of (1) pure dipole, (2) pure quadrupole, and (3) dipolc 
plus quadrupole fields. The quadrupole magnetic moment D is taken to be parallel 
to the dipole magnetic moment /x, and both are inclined relative to the spin axis of 
the star $1 at an angle 0. Simulations have shown that in each case the structure of 
the funnel streams and associated hot spots on the surface of the star have specific 
features connected with the magnetic field configuration. In the pure dipole case mat- 
ter accretes in two funnel streams which form two arch-like spots near the magnetic 
poles. In the case of a pure quadrupole field, most of the matter flows through the 
quadrupole "belt" forming a ring-shaped hot region on the magnetic equator. In the 
case of a dipole plus quadrupole field, magnetic flux in the northern magnetic hemi- 
sphere is larger than that in the southern, and the quadrupole belt and the ring are 
displaced to the south. The stronger the quadrupole, the closer the ring is to the mag- 
netic equator. At sufficiently large 0, matter also flows to the south pole, forming a hot 
spot near the pole. The light curves have a variety of different features which makes 
it difficult to derive the magnetic field configuration from the light curves. There are 
specific features which are different in cases of dipole and quadrupole dominated mag- 
netic field: (1) Angular momentum flow between the star and disc is more efficient in 
the case of the dipole field; (2) Hot spots are hotter and brighter in case of the dipole 
field because the matter accelerates over a longer distance compared with the flow in 
a quadrupole case. 

Key words: accretion, accretion discs - magnetic fields - MHD - stars: magnetic 
fields. 



1 INTRODUCTION 

Many accreting stars have a strong magnetic field which dis- 
rupts the accretion disc and channels accreting matter to the 
star's surface. For such stars many spectral and photometric 
properties are determined by the magnetic field. Examples 
include young solar-type stars (Classical T Tauri stars, here- 
after CTTS, e.g., Hartmann 1994), neutron stars in binary 
systems (X-ray pulsars and millisecond pulsars, e.g. Ghosh 
& Lamb 1978; Chakrabarty et al. 2003), white dwarfs in the 
binary systems (cataclysmic variables, e.g. Wickramasinghe 
et al. 1991; Warner 1995, 2000), and also brown dwarfs (e.g., 
Scholz & Ray 2006). The properties of such accretion will 
depend on the structure of the magnetic field of the star. 
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In early models it was assumed that the star's in- 
trinsic magnetic field was a pure dipole (Ghosh & Lamb 
1979a,b; Camenzind 1990; Konigl 1991). Recently, 2D and 
full 3D MHD numerical simulations of magnetospheric ac- 
cretion were performed which confirmed many predicted 
features and revealed many new details of accretion to a 
star with a dipole field (Romanova et al. 2002, 2003, 2004; 
Long et al. 2005; Kulkarni & Romanova 2005). Different 
aspects of the disc-magnetosphere interaction were also in- 
vestigated numerically (Hayashi et al. 1996; Hirose et al. 
1997; Miller, & Stone 1997; Goodson, Winglee & Bohm 
1997, 1999; von Rekowski & Brandenburg 2004, 2006) and 
theoretically (e.g., Lovelace, Romanova & Bisnovatyi-Kogan 
1995; Uzdensky, Konigl & Litwin 2002; Matt & Pudritz 
2004, 2005). 

However, the actual configuration of the magnetic field 
of strongly magnetized stars may depart from the dipole 
one. For example, Safier (1998) presented a number of ar- 
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guments pointing to a non-dipolar magnetic fields in some 
CTTSs. The Zeeman measurements of the magnetic field of 
a number of CTTSs based on the photospheric lines show 
that the magnetic field at the surface of CTTSs is strong 
(1-3 kGs) but not ordered, which means that close to the 
star the magnetic field is non-dipole (e.g., Johns-Krull et al. 
1999; Johns-Krull, Valenti, & Koresko 1999; Johns-Krull & 
Gafford 2002; Smirnov et al. 2003). 

The measurements of the magnetic fields of rapidly 
rotating low-mass stars with the Zeeman-Doppler imaging 
technique have shown that in a number of stars the magnetic 
field has a complicated multipolar topology close to the star 
(Donati & Cameron 1997; Donati et al. 1999; Jardine et al. 
2002). If the multipolar component dominates in the disc- 
accreting binary systems, then the matter flow to the star 
and the hot spots will be different from those in the case of 
a pure dipole field (e.g., Jardine et al. 2006). 

The dipole magnetic field may of course dominate in 
some stars, and it will dominate at larger distances from 
stars with complex surface magnetic fields. There are obser- 
vational signs that the dipole component possibly disrupts 
the disc. In many cases there is evidence that the dipole field 
dominates at all distances, giving periodic light curves in X- 
ray binaries, intermediate polars, and some CTTSs. In ad- 
dition, direct measurements of the magnetic field using po- 
larization of some spectral lines shows evidence of the pure 
dipole field in a number of the CTTSs (Valenti & Johns- 
Krull 2004; Symington et al. 2005). An remarkable result 
was recently obtained by Donati et al. (2006) who recovered 
the magnetic field of the convective low- mass star v374Peg 
and obtained a pure dipole configuration. Thus the magnetic 
field configuration is likely to be different in stars with vary- 
ing levels of importance of the dipole component. Bouvier 
et al. (2006) argued that there are many signs of magneto- 
spheric accretion in CTTSs and that the field is probably 
dipolar at larger distances from the star, but may have a 
strong multipolar component close to the star. 

Recently, the first theoretical and numerical research 
was done on accretion to a star with a non-dipole magnetic 
field. Jardine et al. (2006) investigated the possible paths of 
the accreting matter in the case of a multipolar field derived 
from observations. Donati et al. (2006) developed a simpli- 
fied stationary model for such accretion, von Rekowski and 
Brandenburg (2006) performed axisymmetric simulations of 
the disc-magnetosphere interaction in case when magnetic 
field is generated by the dynamo processes both, in the star 
and in the disc. They obtained a time-variable magnetic 
field of the star with a complicated multipolar configura- 
tion which shows that the dynamo may be responsible for a 
complex magnetic field structure. 

In this paper we show the first results of full 3D MHD 
simulations of accretion to a star with a non-dipole mag- 
netic field frozen to the stellar surface. As a first step, we 
take a combination of the dipole and quadrupole fields with 
aligned axes and investigate properties of the funnel streams 
and hot spots at the surface of the star. More complicated 
structures of the field will be considered in future simula- 
tions. Analytical analysis of accretion to a non-rotating star 
with a pure quadrupole magnetic field was done by Lipunov 
(1978). Combination of dipole and multipole fields was dis- 
cussed by Lovelace et al. (2005) . However, no numerical sim- 
ulations of such accretion have been performed so far. 



In §2 we describe our numerical model and the magnetic 
field configurations used in the simulations. We describe our 
simulation results in §3, §4 and study the dependence on 
different parameters in §5 and give a summary of our work 
in §6. 



2 NUMERICAL MODEL AND INITIAL 
MAGNETIC FIELDS 

The numerical model used in this paper has been described 
in a number of previous papers (e.g. Koldoba et al. 2002; 
Romanova et al. 2002, 2004; Ustyugova et al. 2006). We de- 
scribe briefly the main aspects of the numerical model in the 
Appendix. In this section, we describe the initial conditions, 
boundary conditions, and new configurations of the dipole 
plus quadrupole field used in this paper. 



2.1 Initial Conditions 

The region considered consists of the star located in the cen- 
ter of coordinate system, a dense disc located in the equato- 
rial plane and a low-density corona which occupies the rest 
of the simulation region. Initially, the disc and corona are 
in rotational hydrodynamic equilibrium. That is, the sum of 
the gravitational, centrifugal, and pressure gradient forces 
is zero at each point of the simulation region. The initial 
magnetic field is a combination of dipole and a quadrupole 
field components which are force-free at t = 0. The ini- 
tial rotational velocity in the disc is close to Keplerian (i.e., 
the small pressure component is taken into account). The 
corona at different cylindrical radii r rotates with angular 
velocities corresponding to Keplerian velocity of the disc at 
this distance r. This initial rotation is assumed so as to 
avoid a strong initial discontinuity of the magnetic field at 
the boundary between the disc and corona. The initial rota- 
tion of the disc with a non-rotating corona leads to a strong 
magnetic forces at the disc-magnetosphere boundary and to 
a strong initial torque of the disc. The distribution of den- 
sity and pressure in the disc and corona and the complete 
description of these initial conditions is given in Romanova 
et al. (2002) and Ustyugova et al. (2006). 

The initial accretion disc extends inward to an inner 
radius r& and has a temperature Td which is much less than 
the corona temperature Td = 0.01T C . The density of the 
disc is 100 times the density of the corona, p c — O.Olpd. The 
values T c , T d , p c , Pd are determined at the boundary between 
the disc and corona near the inner radius of the disc. 

After the beginning of the rotation, the magnetic field 
lines start to twist and the exact force-balance is disturbed 
because the magnetic forces begin to act. However, these 
forces are not strong enough to disturb the disc significantly. 
The accretion rate increases only slightly as a result of the 
magnetic braking associated with this initial twist. These 
initial conditions allow an investigation of accretion to a 
star with a dipole magnetic field for conditions where the 
disc matter accretes inward very slowly on the viscous time- 
scale. The a— viscosity incorporated in our code permits 
the regulation of the accretion rate. Earlier, we used these 
initial conditions to investigate magnetospheric flow in case 
of a pure dipole intrinsic field (Romanova et al. 2002, 2004). 
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2.2 Boundary Conditions 

At the inner boundary (R = R*), "free" boundary con- 
ditions are applied for the density, pressure, entropy, ve- 
locity and </>— component of the magnetic field, dp/dR = 
0, d P /dR = 0, dS/dR = 0, <9(v — fix R,)r/8R = 0, 
d(RB r t > )/8R = 0. The magnetic field is frozen to the sur- 
face of the star, so that the magnetic flux ty(R,9), at 
the inner boundary is derived from the frozen-in condition 
/dt + v p ■ = 0. The poloidal components Br and Bg 
are derived from the magnetic flux function ^(R, 9). 

At the outer boundary, free boundary conditions are 
taken for all variables. Matter does not inflow through the 
outer boundary. The investigated numerical region is large, 
~ 457?*, so that the initial reservoir of matter in the disc 
is large and sufficient for the performed simulations. Matter 
flows slowly inward from external regions of the disc. During 
the simulation times studied here only a very small fraction 
of the total disc matter accretes to the star. Test simulations 
were done where matter was allowed to accrete inward from 
an external boundary did not change the result. 



2.3 Reference Units 





Dimensionless reference values are used in our simulations, 
which allows our models to apply to a range of different sys- 
tems. We let Ro denote a reference distance scale, which is 
equal to 7iU/0.35, where R* is the radius of the star. The 
subscript "0" denotes reference values. We take the refer- 
ence values for the velocity, angular speed and timescale to 
be v = (GM/R ) 1/2 , Q = vo/Ro and P Q = 2nR /v re- 
spectively. The reference magnetic field Bq = B t (R*/Ro) 3 , 
where B, is magnetic field at the surface of the star. The 
reference dipole magnetic moment po — BoRo, and the refer- 
ence quadrupole moment Do = BqRo- The reference density, 
pressure and mass accretion rate are po = Bq/vq, po = PoVq 
and Mo = PoVoR 2 respectively. The dimensionless variables 
are R = R/Rq, v = v/vq, time t = t/Pa etc. In the following 
discussion these dimensionless units are used but the tildes 
are implicit. The grid is Nr = 75, N = 31, which gives our 
external boundary at R ma x ~ 16. 

For a typical CCTS, we take the stellar mass M* = 
0.8Mq, its radius R* = 1.8Rq and the magnetic field at the 
surface of the star B, — lkG, and thus obtain the reference 
units discussed above. Therefore, the case of p, = 0.5 and 
D — 0.5 investigated in the following sections corresponds 
to the magnetic field with the dipole component B*d = 500G 
and quadrupole component B tq — 500G. 



Figure 1. Initial configurations of the magnetic fields used in 
simulations, (a) Pure dipole field with p ^ 0, D = 0; (b) Pure 
quadrupole field D ^ 0; (c) Dipole plus quadrupole field, p ^ 0, 
D 5^ 0. The inner region of the accretion disc is also seen in the 
figures. 



The first term (p^ = because there are no monopoles. 
The second term is the dipole component, = pcosO/r 2 , 



where p is the dipole moment, 
quadrupole component, 



The third term is the 
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where x a , x a are the components of r, D a p is the mag- 
netic quadrupole moment tensor, and the sum over re- 
peated indices is implied. Because D aa = 0, we have, 
^(2) _ D a/ 3n a n{3/2r 3 , where n a — x a /r, np — xp/r. 

Considering the axisymmetric case where Du = D22 = 
—D33/2, we call D = D33 the quadrupole moment and refer 
to the axis of symmetry as the "direction" of the quadrupole 
moment. In spherical coordinates (r,9,<f>), we have, 



2.4 Magnetic Fields of the Star 

In the absence of plasma outside the star we would have 
V x B = 0, so that the intrinsic magnetic field of the star 
can be written as B = — Vip, where the scalar potential 
of the magnetic field is <£>(r) = Em a /|r — r |, and m a is an 
analogy of the magnetic "charge" , r and r a are the positions 
of the observer and the magnetic "charges" respectively. The 
scalar potential can be represented as a multipole expansion 
in powers of 1/r (e.g. Landau & Lifshitz 1983), 
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where P2 is a Legendre polynomial. Then we can get the 
quadrupole magnetic field B and the flux function 
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In our simulations we use a combination of dipole and 
quadrupole magnetic fields: 



B r 
Be 



2/icos6> 3D 2 

— — + i^( 3cos . 

u sin # 3D . . 
-^3— + ^4 sin 61 cos 61 , 



(6) 
(7) 
(8) 



usin 6* 3D . 2 „ 
* = h — r sin cos . 

The star has fixed dipole and/or quadrupole magnetic 
fields with moments /1 and D respectively. For simplicity, 
they are in the same direction and are inclined with respect 
to the rotational axis O at an angle 6. We investigated ac- 
cretion at angles O = 0°, 6 = 30°, and = 60°. Spin axes 
of the star and the disc coincide. 

We vary the value of the quadrupole moment from a 
pure dipole case, D = 0, up to a pure quadrupole case, 
H = 0. Figure 1 shows the main initial configurations of 
the magnetic field considered in the paper: (a) pure dipole 
configuration, D = 0; (b) pure quadrupole case p = 0; and 
(c) the case when the dipole and quadrupole fields are of 
the same order. One can see that when both dipole and 
quadrupole fields are present, the field is not symmetric rel- 
ative to the equatorial plane: the magnetic flux ^ is larger in 
the northern hemisphere where the dipole field is added to 
the quadrupole field, and smaller in the southern hemisphere 
where the dipole field is subtracted from the quadrupole 
field. Simulations have shown that this feature of the dipole 
plus quadrupole field leads to results different from the pure 
dipole case. 



3 RESULTS OF SIMULATIONS 

3.1 Accretion to a Star with Pure Dipole and 
Pure Quadrupole Magnetic Fields 

First, for reference we show the results of simulations of 
accretion to a star with a pure dipole field with /i = 0.5 
inclined relative to the spin axis at an angle O = 30° . Figure 
2a (left panel) shows the density distribution in the x — z 
plane (this is also the \i — CI plane). One can see that most 
of the matter flows through symmetric funnel streams to the 
nearest magnetic pole. A smaller amount of matter flows to 
the opposite pole. The right panel shows the cross- section 
in the y — z plane, where similar but weaker streams arc 
observed. In fact at sufficiently low density levels the matter 
blankets the magnetosphere completely (Romanova et al. 
2003) so that at these levels any cross-section will show the 
funnel streams. 

Figure 2b shows the results of accretion to a pure 
quadrupole magnetic configuration with /i — 0, and D = 0.5, 
which is inclined at an angle 6 = 30° . The left panel shows 
that in the x — z cross-section, matter is slightly lifted above 
the equatorial plane and flows through the quadrupole belt 
to the equatorial region of the star. The right panel shows 
that in the y— z plane, matter flows straight to the equatorial 
belt without being lifted. 

Figure 3a shows the three-dimensional distribution of 
matter flow around a star with a pure quadrupole field, 
where one of density levels is shown in green color. The 
left projection demonstrates that the matter flows in a thin 
sheet approximately in the plane of the magnetic equator. 



The right panel shows that this flow is wide and is not ax- 
isymmetric. Thus, in the pure quadrupole case matter flows 
to the star in a thin, wide sheet in the region of the magnetic 
equator. 



3.2 Accretion to a Star with a Dipole plus 
Quadrupole Magnetic Field 

We calculated the case of the dipole plus quadrupole mag- 
netic field at three misalignment angles, O = 0°, = 30°, 
and O = 60°. In all cases we took /i — 0.5 and D = 0.5. 
This magnetic configuration is not symmetric relative to the 
equatorial plane. In the northern hemisphere the magnetic 
flux of the dipole is added to that of the quadrupole, so 
that magnetic flux is larger than the pure dipole one. In 
the southern hemisphere, the magnetic flux of the dipole is 
subtracted from that of the quadrupole, so that the flux is 
smaller than the dipole flux. 

Figure 4a shows the x — z (left panels) and y — z (right 
panels) cross-sections of accretion to a star with O = 0°. 
One can see that the matter flows axisymmetrically to the 
quadrupole equatorial belt, which is now located in the 
southern hemisphere. 

Figure 4b shows the case O = 30°. Here, part of the 
matter flows to the south polar region. Figure 3b shows the 
3D view of such a flow in the northern hemisphere. One can 
see that matter flows through the quadrupole belt, but the 
flow is not symmetric relative to the z— axis. The right panel 
shows the flow as seen from the southern hemisphere, where 
part of the matter flows to the south pole region. 

Figure 4c shows that at O = 60° , the flow is even more 
complicated: some matter flows to the north pole, some mat- 
ter flows to the south pole, and some matter flows through 
the belt. 

So, one can expect that in the case of O ~ 0°, accre- 
tion to a quadrupole or a combined dipole-quadrupole field 
will go through a quadrupole belt and will form a ring-like 
spot at the surface of the star. However, at large misalign- 
ment angles O the structure of the accreting streams is more 
complicated. 



3.3 Disc-Magnetosphere Interaction and 
Magnetospheric Radius 

Here we chose one of typical cases with mixed dipole and 
quadrupole fields, fx = 0.5, D = 0.5, and O = 30°, and 
consider it in greater detail. First, we estimate the modified 
plasma (3 parameter, 



p + pv 



H B 2 /Stv B 2 /8tv' 

where we took into account that the ram pressure of matter 
in the disc is much larger than the thermal pressure. The 
surface /3 — 1 separates the regions of magnetically dom- 
inated and matter dominated plasma. This surface crosses 
the equatorial plane at the magnetospheric radius r m , which 
is r m ~ 1.14. One can see from Figure 5 that this is the ra- 
dius where the disc is disrupted and matter goes to the fun- 
nel streams. Both dipole and quadrupole components con- 
tribute to the magnetosphere and to the disruption of the 
disc. 
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Figure 2. Disc accretion to a star with (a) pure dipole and (b) pure quadrupole magnetic configurations at = 30° and at time t = 5. 
The left panels show projections in the x — z (fJ.Q) plane, and the right panels show projections in the y — z plane. The background shows 
density contours, and the red lines show the outermost closed field lines. 





Figure 3. Three-dimensional picture of matter flow around the star for the pure quadrupole case (top panels -a ) and for a dipole 
plus quadrupole case at fj, = 0.5, D = 0.5 (bottom panels -b ). In both cases O = 30°. The background shows one of the density levels, 
p = 0.2. Red lines are the outermost closed field lines. In (b), the dashed lines mean the south part of the rotation axis and magnetic 
axis. 



For r < r m , the magnetic energy-density dominates, 
and the field lines are closed. However, at larger radii the 
field lines are dragged by the disc and also innate and open 
in the corona. Figure 5 shows the outermost closed field lines 
(red lines) and also some open field lines (yellow lines) in the 
corona. 

Opening and inflation of the field lines is somewhat 
different between the dipole-dominated and quadrupole- 
dominated fields. Figure 6 shows the large-scale structure of 



the magnetic field in the dipole case at /i = 0.5, D — (left 
panel) and in the quadrupole case at D — 0.5, fi — (right 
panel). One can see that in both cases there is a closed mag- 
netosphere close to the star and a set of inflated field lines at 
larger distance. We should note that in our 3D simulations, 
inflation is somewhat suppressed by relatively dense corona 
and in reality it may be more efficient. From the other side, 
these plots are shown for t — 8 and at longer time inflation 
of the external field lines may progress further. In spite of 
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Figure 4. Disc accretion to a star with a dipolc plus quadrupole magnetic field fi = 0.5, D = 0.5 for different misalignment angles: (a) 
= 0°; (b) = 30°; (c) = 60° at time t = 5. The left panels show projections in the x — z (fJ.£l) plane, and the right panels show 
projections in the y — z plane. The background shows density contours, and the red lines show the outermost closed field lines. 




Figure 5. The detailed structure of disc accretion to a star with a dipole plus quadrupole field, with fi = 0.5, D = 0.5, = 30°. In 
addition to the closed field lines (red lines), a number of sample open field lines (yellow lines) are shown. 



that, the typical pattern of the inflated field lines is differ- 
ent for the dipole and for the quadrupole. In both cases the 
open field lines participate in the spin-down of the star (see 
Section 5). 



4 HOT SPOTS AND LIGHT CURVES 

The inflowing matter of the funnel streams falls onto the star 
and forms hot spots on the stellar surface. The shape of the 
hot spots reflects the shape of the funnel streams. Earlier 3D 
simulations of accretion to a star with a misaligned dipole 
magnetic field have shown that hot spots are always simi- 
lar in the northern and southern hemispheres (Romanova 
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Figure 6. The large-scale structure of the magnetic field (yellow lines) in the cases of pure dipole at p = 0.5 (left panel) and pure 
quadrupole D = 0.5 (right panel) at = 0. The background shows the density distribution where density varies from p = 2 (red) to 
p = 0.01 (blue). 




Figure 7. The hot spots at t=5 shown at different projections: from the top to the bottom, each column shows views from the north 
and south magnetic poles, the x — z projection and the p — y projection respectively. The background shows density contours. Red color 
corresponds to the maximum density (p = 2.5) and dark blue corresponds to the minimum density p = 0.01. (a) pure dipole case, p = 0.5, 
D = at = 30°; (b) pure quadrupole case, p = 0, D = 0.5 at = 0°. Columns (c), (d) and (e) show spots for the p = D = 0.5 case 
for = 0° , = 30° and = 60° respectively. 



et al. 2004). Our present simulations with a dipole plus 
quadrupole magnetic fields have shown that hot spots are 
always different in the northern and southern hemispheres. 

Figure 7 shows the hot spots for different magnetic con- 
figurations at t — 5. From the top to the bottom, each col- 
umn shows view from the north and south magnetic poles, 
the x — z projection and the /i — y projection respectively. 

Column (a) represents the results for a pure dipole mag- 



netic field with jj, — 0.5 and misalignment angle O = 30°. 
One can see that there are two hot spots near the magnetic 
poles which are similar to each other. Column (b) shows the 
results for a pure quadrupole configuration at O = 30° . The 
inflowing matter forms a ring near the magnetic equatorial 
plane in the region of the quadrupole belt. The ring has two 
density enhancements from two opposite sides because the 
flow is not axisymmetric (See also Figure 3a). Columns (c), 
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(d) and (e) show the results for a dipole plus quadrupole 
field, fx — 0.5, D — 0.5, for misalignment angles O = 0, 
= 30° and = 60° respectively. In the case of = 0°, 
(column c), the hot spot has the shape of the ring located 
below the equatorial plane. The ring is below the equator 
because the magnetic field in the northern hemisphere is 
much stronger than that in the southern hemisphere and 
the quadrupole belt is displaced to the south. The ring is 
symmetric relative to ft axis because = 0°. There are no 
hot spots near the poles (see also Figure 4a). Column (d) 
shows the case when = 30°. One can see that the mag- 
netic flux in northern hemisphere is still strong enough to 
stop the flow to the north magnetic pole, so there is no hot 
spot. However, some matter flows to the south magnetic pole 
along closed field lines of the southern part of the magneto- 
sphere. Column (e) shows the results for O = 60°. Now the 
magnetic axis is tilted more towards the disc, which leads 
to matter flow to both the north and south magnetic poles. 
We can still see the belt in the southern hemisphere. 

Consider now the light curves from hot spots on the 
star for different magnetic configurations. We assume that 
the kinetic energy of the inflowing matter is converted into 
black-body radiation and is radiated isotropically (e.g. Ro- 
manova et al. 2004). The variation in the light curve is as- 
sociated with the position and shape of the hot spots. This 
means that the evolution of the hot spots can affect the light 
curves. After several rotations of the star, the shape and po- 
sition of hot spots do not vary significantly. Therefore, as a 
first approximation we calculate the light curves by fixing 
the hot spots at t = 5 and rotating the star. 

The observed intensity of radiation in the direction k 
is J = J/(R, k)dS, where /(R, k) is the intensity of the 
radiation from a unit area into the solid angle element dQ, 
in the direction k, dS is an element of the surface area of 
the star, R is the radius vector. It can be calculated from 
our simulations (see details in Romanova et al. 2004). 

Figure 8 shows the light curves for a pure dipole, pure 
quadrupole and a dipole plus quadrupole case at = 30° at 
different inclination angles from i = 15° to 90°. The inclina- 
tion angles are measured relative to the star's rotation axis. 
One can see that for the pure dipole case (top panel) the ob- 
server will see only one peak per period, excluding the case of 
very high inclination angles, i ^ 70°. In the pure quadrupole 
case (middle panel), the observer will see two peaks per pe- 
riod starting from smaller inclination angle, i ^ 40°. In the 
case of the mixed dipole plus quadrupole field, there are also 
two peaks per period at i ^ 40° . Also, when the inclination 
angle is changed, the change in amplitude is stronger in this 
case. If the inclination angle of the star is known, one can 
derive the magnetic configurations by varying O. However, 
there may be several possible configurations of the field and 
which give similar light curves. 



5 DEPENDENCE ON THE STRENGTH OF 
THE QUADRUPOLE COMPONENT 

Here, we consider the dependence of different features of 
the accretion flow and hot spots on the strength of the 
quadrupole component. We fixed the dipole component at 
ft — 0.5 and varied the quadrupole component between 
D = and D = 2.5. First, we calculated the mass accre- 
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Figure 8. Light curves for different magnetic field configurations 
at misalignment angle G = 30° for different inclination angles i. 
The inclination angle is measured relative to the star's rotation 
axis. Solid, dash, dash-dot and dotted lines represent the results 
for i = 15°, i = 30°, i = 60°, i = 90° respectively, (a) Pure dipole 
case fi = 0.5, D = 0; (b) Pure quadrupole case fi = 0, D = 0.5; 
(c) dipole plus quadrupole case at fi = 0.5, D = 0.5. 



tion rate M = — J dS ■ pv p along the surface of the star S 
and observed that accretion rate varies by different way but 
on average it does not depend on the value of the quadrupole 
component (see Figure 9a). This is an expected result, be- 
cause the accretion rate should depend on the viscous prop- 
erties of the disc but not on magnetic configurations of the 
star. All simulations were done at the same parameters of 
the disc and fixed value of a— viscosity, a = 0.04. We should 
note that initially, during the first 1-2 rotations accretion 
rate is very small. This is the time when matter moved to- 
wards the disc and reached the surface of the star through 
the funnel streams. Later, the structure of the spots estab- 
lishes and varies only slightly with time. At present, we are 
able to simulate accretion to a star with a quadrupole field 
during 5 — 50 rotations. We show the plot in Figure 9 up 
to t = 5 because the spots do not change appreciably for 
longer times. 

We also calculated the spinning torque to the star 
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Figure 10. Figure shows streamlines of angular momentum flux carried by the magnetic field lines fs = rB^Bp/An in the cases of the 
pure dipole field at /i = 0.5 (left panel) and pure quadrupole field D = 0.5 (right panel) at G = 0. The background shows the angular 
velocity distribution where angular velocity varies from Q = 0.4 (red) to O = 0.01 (blue). 
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Figure 9. Mass accretion rate M (a) and angular momentum 
fluxes L f (b) for stars with different strength of quadrupole mo- 
ment which changes from D = to D = 2.5, while fi = 0.5, 
= 30° and other parameters are fixed. 



from the disc-magnetosphere interaction. The main spinning 
torque comes from the magnetic field lines connecting a star 
with the disc and corona Lf = J dS ■ rB$B v j (Atk) (see also 
Figures 5 and 6) . A star rotates relatively slowly (co-rotation 
radius is r cr = 2) and magnetospheric radius is smaller than 
co-rotation radius (r m ~ 1) so that the inner regions of the 
disc rotates faster than the star and accreting matter spins 
the star up through the closed field lines connecting a star 
with the disc. 

There is also a negative torque from the open field lines 
which start in the polar regions and connect a star with 
slowly rotating corona. This torque leads to spinning-down, 
but this torque is smaller compared to the spin-up torque. 



Figure 9b (solid line) shows the torque in the case of the pure 
dipole star. When we added small quadrupole component, 
then torque stayed approximately the same. However, when 
we added the largest quadrupole component (D = 2 and 2.5) 
then the positive torque decreased significantly. Analysis of 
angular momentum fluxes in the pure dipole and quadrupole 
cases have shown that in the dipole case the main positive 
torque comes through the inner closed field lines which con- 
nect the disc and the star. Disc rotates faster than the star 
and angular momentum flows to the star. In the quadrupole 
case disc matter flows through the equatorial belt and there 
is less "communication" between the star and the disc. This 
is why in quadrupole case the positive torque is smaller. This 
explains the fact that at very large quadrupole component 
the spinning torque is smaller: the field structure becomes 
very similar to a pure quadrupole topology and the disc ac- 
cretes through a quadrupole belt. 

We investigated the angular momentum flow in greater 
detail in the cases of pure dipole field (/i = 0.5, D = 0) 
and pure quadrupole field (D = 0.5, fi — 0). Figure 10 
shows the distribution of the angular velocity (background) 
and streamlines of the angular momentum fluxes carried 
by the magnetic field /b = rB^Bp/Air. One can see that 
in the dipole case, angular momentum flows in from the 
disk to the star (the star spins-up), and flows out through 
the open polar field lines into the corona. In the case of 
the quadrupole field, angular momentum flows to the star 
from the quadrupolar equatorial belt, where the disc mat- 
ter accretes and is threaded by some closed field lines of 
the quadrupole. This is the region where the inward angular 
momentum flux is the largest. Again, there is angular mo- 
mentum outflow from the star through the open polar field 
lines. In both cases the total angular momentum transferred 
to the star is positive and the star spins-up. We should note 
that the angular momentum carried to the star by matter is 
much smaller than that associated with the field, by a fac- 
tor of 20-100 in the dipole case and 10-20 in the quadrupole 
case. 

In addition we investigated different properties of the 
hot spots in cases of aligned (O = 0°) pure dipole and 
pure quadrupole configurations at equal magnetic moments, 
H — 0.5 and D = 0.5. In the case of an accreting dipole 
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the matter of the funnel stream is lifted above the mag- 
netosphere and then is accelerated towards the surface of 
the star, so that significant part of gravitational energy is 
released thus increasing the temperature of the hot spots 
(see Romanova et al. 2004). In the case of predominantly 
quadrupole field matter flows through the equatorial belt 
without such dramatic falling to the surface of the star. The 
simulations show that the velocity of the flow and the Mach 
number are about two times smaller in the quadrupole case. 
The luminosity of hot spots depends on the density and ve- 
locity. The density is similar in both dipole and quadrupole- 
dominated cases, difference in velocities leads to conclusion 
that the hot spots are expected to be (vd/v q ) 2 times dimmer 
in the quadrupole case. The temperature of the hot spots is 
expected to be smaller in the quadrupole case by the factor 
(vd/v q ) 3 ' 14 ' (about 1.7 times) which follows from equating of 
the kinetic energy flux to the black-body radiation flux. In 
both cases the maximum energy and highest temperature 
are in the middle of the spot. Thus, a star with pure dipole 
case is expected to have hotter and more luminous spots. 



6 SUMMARY 

We investigated for the first time disc accretion to a star 
with a dipole plus quadrupole field in full 3D MHD simu- 
lations. The simulations show that for relatively small mis- 
alignment angles 0, matter flows in a thin wide sheet to the 
quadrupole's belt. This flow forms a ring-like hot spot at the 
surface of the star. The position of the ring coincides with 
magnetic equator in the case of the pure quadrupole field, 
and is displaced into the southern hemisphere in the case of 
the dipole plus quadrupole fields. At large 6, most of matter 
may accrete to one or both magnetic poles. In quadrupole 
case the light curves typically show two peaks per period 
for the majority of inclination angles i. This might be a 
sign of the significant quadrupole component. However in 
the case of pure dipole field two peaks per period are typi- 
cal for the case of high inclination of the system (Romanova 
et al. 2004). If inclination of the systems is known, then 
the light curves may help to suggest the possible role of the 
quadrupole component. 

The simulations show that the torque on the star is 
larger in the case of the dipole field, where angular momen- 
tum is transferred through the closed field lines connecting 
the star and the disc. In the quadrupole case the torque 
is smaller because matter of the disc accretes directly to 
the star through the quadrupole belt, that is, between field 
lines. There is a positive torque associated with the disc- 
magnetosphere interaction, however it is much smaller than 
in pure dipole case. 

In the case of the dominant quadrupole field, the hot 
spots are expected to be cooler with smaller total energy 
released because in the quadrupole belt the funnel flow hits 
the star with smaller velocities compared to the pure dipole 
case. To apply these predictions, one needs to know the ac- 
cretion rate and the value of the magnetic field on the surface 
of the star. 

In reality the star's magnetic field may be even more 
complicated and may include higher order multipoles. In 
such a field multiple funnel streams are expected to form 
(Donati et al. 2006; Jardine et al. 2006). Such multipolar 



fields were recently modeled by von Rekowski and Branden- 
burg (2006) in the axisymmetric time-dependent simulations 
which incorporated the dynamo processes of the field forma- 
tion both in the disc and in the star. In the case of multi- 
polar fields the light curves are expected to be even more 
complicated. The light curves obtained in this paper and 
the light curves obtained earlier for pure dipole field (Ro- 
manova et al. 2004) may help to distinguish the cases where 
the magnetic field is highly ordered (dipole or dipole plus 
quadrupole) from the cases of the multipolar field where the 
light curves are expected to be much less ordered. In the 
future research we plan to investigate the case where the 
dipole and quadrupole are misaligned and also the case of 
accretion to a star with higher order multipolar magnetic 
field. 



ACKNOWLEDGMENTS 

This research was conducted using partly the resources of 
the Cornell Theory Center, which receives funding from Cor- 
nell University, New York State, federal agencies, founda- 
tions and corporate partners, and partly using the NASA 
High End Computing Program computing systems, specif- 
ically the Columbia supercomputer. The authors thank 
A.V. Koldoba and G.V. Ustyugova for the earlier devel- 
oping of codes and A. K. Kulkarni for helpful discussions. 
This work was supported in part by NASA grants NAG 
5-13060, NNG05GL49G and by NSF grants AST-0307817, 
AST-0507760. 



APPENDIX 

We perform full 3D MHD numerical simulations. Equations 
are written in a reference frame rotating with the star, with 
the z— axis aligned with the star's rotation axis, 



|+V.(p«) = , 



djpv) 
dt 

d(pS) 

dt 
dB 
dt 



+ V • T = pg + 2pv x SI - pSl x (SI x R) 



+ V ■ (pSv) = 



V x (v x B) = , 



(9) 



(10) 



(11) 



(12) 



where v and B are vectors of velocity and magnetic field in 
the three-dimensional space, p is density, St is the angular 
velocity of rotation of the star, S is the entropy per gram, 
g is the gravitational acceleration, T is the stress tensor. 
Using the reference frame rotating with the star we eliminate 
from equation large terms corresponding to the magnetic 
forces induced by rotation of the strong field associated with 
connected with rotation 

The problem of accretion to a star with a dipole and 
quadrupole magnetic fields is complicated because of high 
gradients of the magnetic field, B d ~ 1/R 3 for the dipole 
and B q ~ 1/R 4 for the quadrupole components. To avoid 
this difficulty and other difficulties connected with strong 
magnetic field near the star, the magnetic field was decom- 
posed into the "main" component, which is fixed and con- 
sists of the dipole and quadrupole parts: B = Bd + B q , and 
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variable component B± which is induced by currents in the 
simulation region and is calculated in equations. This split- 
ting helped to decrease the magnetic force (Tanaka 1994; 
Powell et al. 1999). 

A special "cubed" sphere grid was developed for solu- 
tion of such a problem which has advantages of the spherical 
and Cartesian coordinate systems (Koldoba, et al. 2002). A 
grid consists of a sequence of co-centric spheres (cubes in- 
flated to the sphere), which are placed at different radii R 
from the center of coordinates (where a star is located). At 
each inflated cubed sphere there are six sectors correspond- 
ing to the sides of the cube, and curvilinear Cartesian coor- 
dinated are introduced in each of six sides with TV 2 cells. The 
concentric "spheres" are inhomogeneously distributed such 
that the grid is very fine closer to the star and it is much 
coarser far from the star. Each grid cell has approximately 
the same sides in all three directions and this corresponds 
to refine grid near the star. Thus, the whole simulation re- 
gion consists of six blocks with Nr x TV 2 cells. In the current 
simulations we took the grid 75 x 31 2 in each of six blocks. 
Other grids were investigated for comparison. The coarser 
grids give satisfactory results in the case of the pure dipole 
field, however in the case of the quadrupole field, the code 
requires finer grid, because of the higher magnetic field gra- 
dients. The 3D MHD code was parallelized, so that the jobs 
can be run on N = 6 x n processors. At TV > 60 the commu- 
nication time between R— regions becomes significant, but 
in the future in the case of finer grid, we can use even larger 
number of processors. 

To solve the above partial differential equations, we used 
a high order Godunov-type numerical code earlier developed 
in our group (Koldoba et al. 2002). The numerical scheme 
is similar to one described by Powell et al. (1999). The 
Godunov-type schemes have been recently well described in 
the literature (see, e.g. book by Toro 1999). In addition we 
incorporated viscosity to the numerical scheme to control 
the rate of matter flow to the star. The main goal of vis- 
cosity is to bring matter from the far regions of the disc to 
the star, so that we used a simplified a— prescription with 
a = 0.04 in all simulation runs with smaller/larger values 
for testing. 
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